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Abstract 


Using exact formulae for the scattering data of the Benjamin-Ono equation valid for general 
rational potentials recently obtained in llT9l . we rigorously analyze the scattering data in the 
small-dispersion limit. In particular, we deduce precise asymptotic formulae for the reflection 
coefficient, the location of the eigenvalues and their density, and the asymptotic dependence of 
the phase constant (associated with each eigenvalue) on the eigenvalue itself. Our results give 
direct confirmation of conjectures in the literature that have been partly justified by means of 
inverse scattering, and they also provide new details not previously reported in the literature. 
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1. Introduction 

The Benjamin-Ono (BO) equation 



—r = 0, —oo < X < oo, f > 0 

dx^ 


( 1 ) 


describes the weakly nonlinear evolution of one-dimensional internal gravity waves in a stratified 
fluid |ITl|4l|2T], where u corresponds to the wave profile, e > 0 is a measure of the effect of 
dispersion, and the operator “K denotes the Hilbert transform defined by the Cauchy principal 
value integral 



( 2 ) 


The BO equation ([T]i has been used to model internal waves in deep water Q, the atmospheric 
roll cloud wave-train known as the morning glory ll2^ . and nonlinear Rossby waves in shear 
flow II 22 I . to name a few examples. Beyond fluid dynamics, equation Q has been observed to 
model the spectral dynamics of incoherent shocks in nonlinear optics lH. 

In the small-dispersion limit (e —> 0), numerical experiments indicate that, for t sufficiently 
small (independent of e), the solution of the BO equation with smooth e-independent initial data 
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Mo is well approximated by the solution of the inviscid Burgers equation (Hopf equation) ob¬ 
tained by setting e = 0 in ([1]). At the advent of a shock in the dispersionless model, the solution 
of (0 is then regularized with the formation of a dispersive shock wave (DSW). This DSW has 
the well-known structure of an (9(l)-amplitude modulated periodic traveling wave with wave¬ 
length 0(e). In the DSW region the solution of the BO equation may be formally approximated 
using Whitham modulation theory. Surprisingly, and unlike the case of the Korteweg-de Vries 
(KdV) equation, the modulation equations for the BO equation are fully uncoupled Q, consist¬ 
ing of several independent copies of the inviscid Burgers equation. A formalism for matching 
the Whitham modulation approximation for the DSW onto inviscid Burgers solutions in the do¬ 
main exterior to the DSW was developed and applied by Matsuno iniiisi and Jorge, Minzoni, 
and Smyth ||9l to analyze the Cauchy problem for ([T]). Partial confirmation of these results was 
given by Miller and Xu Il20ll . who rigorously computed the weak limit of (essentially, modulo an 
approximation of the scattering data) the solution of the Cauchy problem for (0 for a class of pos¬ 
itive initial data using inverse-scattering transform (1ST) techniques, by developing an analogue 
for the BO equation of a method first invented for KdV by Lax and Levermore II2. Avoiding 
the (a priori unjustified) approximation of the scattering data requires its careful analysis in the 
small-dispersion limit by direct means. 


One application of small-dispersion theory in nonlinear waves involves the investigation of 
universality, for example determining properties of general solutions near the onset of the DSW 
(the point of gradient catastrophe for the inviscid Burgers equation) that are independent of ini¬ 
tial data. It has been postulated by Dubrovin |j6l that the solutions of a wide class of dispersive 
Hamiltonian perturbations of the inviscid Burgers equation are universally modeled by a partic¬ 
ular solution of an integrable fourth order ODE (Painleve I 2 ) near the gradient catastrophe point. 
This result was subsequently proven rigorously for the KdV equation by Claeys and Grava a 
using Riemann-Hilbert analysis in an appropriate double scaling limit. Interestingly, the BO 
equation does not fall within the Dubrovin universality class of equations. However, Maseoro, 
Raimondo, and Antunes lfT4l have recently generalized Dubrovin’s method to the BO and other 
nonlocal equations. Based on this analysis, in IIT4l it is postulated that the solution of the BO 
Cauchy problem should be universally modeled near the point of gradient catastrophe by a par¬ 
ticular solution of a certain nonlocal analogue of Painleve I 2 . Rigorous verification of this formal 
result using the 1ST would again require, among other tools related to asymptotic analysis of non¬ 
local Riemann-Hilbert problems yet to be developed, the accurate asymptotic characterization of 
the scattering data for the BO equation. 


In this paper we use recently obtained exact formulae in to study the asymptotic behavior of 
the scattering data for the BO equation ([T]) in the limit e ^ 0. These formulae hold for rational 
initial data with simple poles, and for convenience we adopt here the further condition of a single 
local extremum (Definition 3.11. The formulae for the scattering data are sufficiently explicit 
to admit rigorous asymptotic analysis by classical methods such as steepest descent (for inte¬ 
grals), thus justifying and extending previously reported formal asymptotic results. We confirm 
directly the celebrated formulae of Matsuno for the asymptotic density of the eigenvalues m 
(Corollary |5.2| l and the magnitude of the reflection coefficient ifT^ (Corollary |4.1| l, the formula 
postulated in Miller and Xu Il20l for the asymptotic values of the phase constants (Theorem |5.2| i, 
and Xu’s conjecture on the phase of the reflection coefficient with some modifications (The- 
orem |4.l| i. We also obtain more precise information about the discrete spectrum than had even 
been conjectured before (Theorem [5T] and Corollary |5.1| l. 
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2. Direct Scattering for Rational Initial Data 

The direct scattering problem for the BO equation can be thought of as the problem of spectral 
analysis of the operator £. := -iedx + C^mqC^ acting on the Hardy subspace H^(W) of L^(IR) 
consisting of functions analytic in the upper-half ;c-plane. Here, the potential mq : K ^ M is the 
initial condition u{x, 0) = uq{x) for Q, and the operator : L^(M) ^ is the self-adjoint 

orthogonal projection onto the Hardy space. As first noted in ifTTI and generalized in IT9l . the 
direct scattering problem for the BO equation can be effectively solved for rational potentials of 
the form 

p 

V 1 ^ P 1 

m(x, 0) = Mo(x) = >-h c.c. satisfying mq 6 L (M), (3) 

^ X- Zn 
p=i ‘ 

where Cp and the poles are points with positive imaginary parts and distinct real parts 

increasing with p. The condition uq e L^(K) is equivalent to Rejci -H ■ ■ ■ H- cp] - 0, in which 
case M := uo(x) dx - 2m(ci H- ■ ■ ■ H- cp) e K. For potentials of the form Q, the problem of 
recovering the scattering data is essentially reduced to the problem of studying the linear system 

A(d)v(/i) = b(d) (4) 

for an unknown vector v(/l) e C^, where A(/l) e and b(/l) 6 have components 

- dz and bJA):^-A dz (5) 

C,„ z - Zp Jc„ 

in which the exponents involve the function 

h{x'. A) := Ax + f(x), (6) 

and / denotes a certain anti-derivative of uq: 

p 

fix) ^ Cp ^Log (/(x - Zp)) -t- y j H- C.C., fix) = mo(x), (7) 

p=i 

where Log(-) is the principal branch (|Im{Log(-))| < tt). The contours Ci,..., Cp lie in the domain 
of analyticity of / and are chosen as follows when T < 0. If ic,„f is not a negative integer, then 
Cm = Uf a contour beginning at ioo to the left of the line Re{z) = Rejzi) and terminating at ico 
to the right of the line Re{z) = Rejzm) and (if m < P) to the left of the line Re{z) = Re{z„,+i). 
If icmf is a negative integer, then C„, = toiZm) is a contour beginning at ioo to the left of the 
line Re{z) = Rejzi) and terminating at Zm- See Figure [T] The matrix A(/l) and vector b(/l) can 
be analytically continued from T < 0 to the domain C \ and while the formulae (|^ remain 
valid, due to the dominant term Ax in the function h, the contours Ci,..., Cp must be suitably 
rotated (continuing / through its branch cuts) to maintain absolute convergence of the integrals 
and allow the continuation. 

The discrete eigenvalues of the BO spectral problem for the rational potential uq are the roots 
of the Evans function 

DiA) det(A(T)), T 6 C \ K+. (8) 

As £. is self-adjoint, the eigenvalues are real numbers. Moreover, since -ied^ is positive semi- 
definite on H^iW), it is easy to see that the eigenvalues A satisfy 


A> - sup{Mo(-*)), 4 an eigenvalue. 

xeR 
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( 9 ) 





Figure 1: The domain of analyticity of the function / omits from C vertical branch cuts emanating from each of the poles 
of Wo in the direction away from the real axis, illustrated here for a rational potential wq with P = 3. Left: contours 
Right: contours 


Associated with each eigenvalue A - Aj < 0 there is a unique normalized (although not in the 
L^(M)-sense) eigenfunction Oy(x) = <I>(x; Aj), where (convergent improper integral) 

p 

^ r y tAAL dz. (10) 

t J~oo Z Zp 

Here, the coefficients ^p(A) are defined for all d e C \ by 

p 

Y^A„,p{A)(l)p{A)^0, ( 11 ) 

p=i 


together with 

p 

Y^(Pp(A)^A. ( 12 ) 

/7=I 

When /I = /ly < 0 is an eigenvalue, A(/l) has rank P - and thus ^{A) - (0i(/l),..., is 

the unique null vector of A(/lj) normalized by the condition ( [T^ . 

The scattering data associated with the rational potential mq are as follows HD. 


Continuous spectrum for T > 0. 

• Reflection coefficient [i{A) for T > 0: defined by 


m- 




r 


-ih(zU)l€ 


V„(d)' 




P=1 


dz 


(13) 


with M 2m{c\ + ■■■ + cp) e M. The coefficients {Vp(d))^^[, assembled in a vector 
v(/l) := (vi(A),..., vp(A))^, satisfy the linear algebra problem 


A^(T)v(/l) = b>(T), T > 0, 


(14) 


where A^(/l) 6 and b^(/l) 6 denote the boundary values taken by A(/l) and b(/l), 
respectively, on from the upper-half d-plane. 
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Discrete spectrum for /I < 0. 

• Negative eigenvalues the zeros of the Evans function D{A) defined by ([^. 

• Phase constants defined by 

X CXD 

(xO/x) - l) dx (15) 

oo 


in terms of the corresponding eigenvalue Aj and normalized eigenfunction <l>j(x) = <l)(x; Aj) 
defined by ( [TOl i. 


Remark 2.1. The number N of eigenvalues Aj and phase constants jj is dependent on the disper¬ 
sion parameter e. In fact, N — N(e) — as e ^ Ofor P fixed; see Section 5.1 


Remark 2.2. The formula ( |15| l is a relation known to hold for general real potentials uq /[T^. 
In l\19\l a different formula for jj was shown to be valid for rational potentials uq. While that 
formula is especially useful for exact calculations (see Section^, the general formula ( |15| l turns 
out to be better suited to small-dispersion asymptotics; see Section\5^ 


3. Rational Klaus-Shaw Potentials 

The analysis of the BO scattering data in the small-dispersion limit is the most straightforward 
in the case that the rational potential mq is of Klaus-Shaw type. 

Definition 3.1. A rational Klaus-Shaw (rKS) potential uq is a rational function of the form 0 
for which there exists Xc e K such that uq is strictly monotone on the intervals (—oo, xf) and 
(Xc, -l-oo) (making Xc the unique real critical point). 

Rational Klaus-Shaw potentials are of one sign, and hence there are distinct classes of positive 
and negative rKS potentials. Colloquially one can describe a rKS potential as having a graph 
consisting of a single “lobe” or “bump,” a property that was found by Klaus and Shaw Col to be 
useful in confining the spectrum of certain non-self-adjoint operators. This property also provides 
an important simplification in various asymptotic calculations involving the WKB approximation 
for differential equations with Klaus-Shaw coefficients, for which there are at most two real 
turning points. The utility of the monotonicity condition is similar in the context of the BO 
scattering problem in the small-dispersion limit. Indeed, for a rKS potential uq there are at most 
two real roots T e K of the equation mo(x) = -A, and hence the corresponding inverse function 
x(A) has two real branches. This implies some useful properties of the exponent function h 
defined in (|^, because the critical points of h satisfy 

h'(x; A) - A-h uo(x) - 0 or uo(x) = -A; (16) 

here the prime denotes a derivative in x. 

Definition 3.2. The bulk S c is the set S - {A : - supjwo) < T < - infjMo))- We say that A e S 
lies in the bulk, while T e K \ S lies outside the bulk, where S denotes the closure ofS. 

Remark 3.1. For a strictly positive (resp., negative) rKS potential uq, S — (— supjwolj 0) ( resp., 
S — (0, — inf{Mo))) so the bulk is a negative (resp., positive) interval abutting the origin. We refer 
to the origin as the “hard edge” of the bulk and to the nonzero endpoint of IB as the “soft edge” 
of the bulk. 
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Definition 3.3 (Critical points for A in the bulk). Let uq be a rKS potential with corresponding 
bulk IB. The functions x± : B —> M represent the two real branches of the inverse of —Uq with 
X-(A) < x+(A)for A e B and are the only real critical points of the function h; see Figure^ The 
remaining 2P — 2 critical points of h form complex-conjugate pairs for A e B and will be denoted 
by Xp, Xp with lm{Xp} > 0 for p — I,..., P — I (not necessarily distinct for all A € B). 



Figure 2: A strictly positive rKS potential uq illustrating the two real branches jr±(/l) of the inverse function of -uq 
defined for suitable A < 0. For a strictly negative rKS potential the picture is similar, but i > 0 is necessary for the 
functions Jc±(/1) to be defined. 


Definition 3.4 (Critical points for A outside the bulk). Let mq be a rKS potential with correspond¬ 
ing bulk B. If A & M.\B, then h has 2P critical points in complex-conjugate pairs denoted by 
Xp, x*p with \m{xp] > Qfor p — 0,... ,P — I (not necessarily distinct for all A e.^\B). 


In a neighborhood of d e S or d 6 K \ S for which all 2P critical points of h are simple, each 
critical point is a locally analytic function = x(A) of d that satishes ([T^. As d exits the bulk 
through the soft edge, the two real critical points x±(A) coalesce and bifurcate into the complex 
plane, where they become renamed as the (nearly real) conjugate pair jco(d) and jco(d)*. As d exits 
the bulk through the hard edge, the two real critical points x±{A) tend to +oo and then become 
finite again as a (large) conjugate pair ;co(d) and xo(d)*. 

Since according to Q and h'(x', d) is a rational f uncti on of x with simple poles at 
such that h'{x', d) —> d as x — > oo, using Dehnition 

h'(x- A) = d'l'^(x; d)'l'+(x; d) for 


3.3 we may explicitly write 
AeB, (17) 


where, assuming the complex critical points Xp = Xp(A) are distinct for A e B, 

np=i ~ Xp(A)^ 


'l'^(x;d) := (x-x-(d)) 

'l'+(x;d):=(x-x+(d)) 


0^=1 {x - Zp) 


( 18 ) 


Remark 3.2. Note that ('k^j is analytic in the upper-half (lower-half) x-plane. For A outside 
the bulk, we may also write ( |17[ ) provided we extend the definitions of'P^ by replacing x_ with xq 
and x+ with Xq in These extended definitions satisfy 'k (x; d)* = ^^^(x*; A) for all A e]S.\B. 
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Differentiating (O and using the definitions ([TSll, we may express d) for A e S as 


ivur (y ("tv 

= = for AeS. (19) 

From equation ( [T9l ) we may immediately observe that 

sign{/!"(x±(/l); d)) =+sign{d) for AeS. (20) 

This fact can also be deduced from equation ( |T6] l and Figure Namely, for a positive rKS 
potential uq, we have u'f^(x-(A)) > 0 and Uq(x+(A)) < 0 for a given d e S, so /r"(x±(d);d) = 
+|Mg(;ic±(d))|. Similarly, for strictly negative uq, we have that h" (x±(A); A) = +|Mg (;ii:±(d)) | for 
AeS. 

We close this section by showing how to evaluate the function h(x; d) at its real critical points 
for d 6 S. 

Lemma 3.1. Let A e S, and set h±(A) h(x±(A); A). Then 


h^{A) = 6»+(d) + M and h^A) = 0-(d), 


( 21 ) 


where 

X c» pA 

uo(y)dy and 0±(A) I x±{rj)dri 

OO v/0 

are well-defined because uq 6 L^(R). 

Proof. Because f{x) is an antiderivative of uq{x) that vanishes as jc —> -oo. 


( 22 ) 


fixAA))^ I uo(y) dy ^-AxAA) - I yu[fy)dy, 

V/— OO %J—OO 


(23) 


by integration by parts. We rewrite the final integral on the right-hand side using the substitution 
Mo(y) = -t], where y - xfirj) for y 6 (-oo, x_(d)]. Hence equation (|2^ becomes 


f 


f(x-(A)) = -Ax-(A) + I x-(7j)dj], 
i.e., /j-(d) = 0_(d). Starting instead from the formula 


X OO ^OO 

uo(y) dy- i 

OO v/a;+(/}) 


/(.x+(d))= uo(y)dy 


uo(y) dy 


(24) 


(25) 


but making the substitution y - x+if) after integrating by parts in the second term then gives the 
desired formula for/!+(d). □ 


4. Continuous Spectrum d > 0 

An asymptotic expansion for the reflection coefficient y6(d), d > 0, can be directly obtained 
from formula ( [T3] l and the linear system ( [T4l i. We first prove that the solution v(d) of the linear 
system ( [T4l i has a limit v‘’(d) as e ^ 0, at which point it only remains to apply the method of 
stationary phase to the integral in ( [T3] l to determine the leading-order behavior of jS(d). 
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Every term in the m* equation of the linear system ( [T4l i is an integral over a common contour 
Cm that is obtained from (appearing in the definition of the system Q) by rotation of the 
vertical tails through n radians. The multi-valued factor g appearing in the integrand 

is analytically continued through its vertical branch cuts as the contours rotate so that when 
T > 0, A^(/l) and b^(/l) are the boundary values taken on from the upper half-plane of the 
corresponding quantities A(/l) and b(/l) analytic for /! e C \ We may then define a related 
linear system in which each row is first taken as a linear combination over Z of the rows of 
and then simplifications accounting for factors arising from continuation of e about its 

singularities made. This new system has the same form as ([T^ but each integral over 

Cm is replaced by an integral of some analytic branch of the same integrand over a new contour 
which we denote by Wm, and it can be written in the form NA^(/l)v(/l) = Nb^(/1) for some 
square matrix N. Our asymptotic analysis of the solution v(/l) of the system ( [T4l l will be effective 
provided the new system is suitable in the sense that (i) each contour Wm passes through exactly 
one critical point of h(z', T) such that Re{-ih{z', d)) is maximized along Wm at that point, (ii) the 
association of critical points to contours W,„ is one-to-one, and (iii) the matrix N is invertible 
(this property guarantees the equivalence of the old and new systems). 

We believe that in general a suitable linear system can be found. We present the following 
proposition as proof of principle in a general case. 


Proposition 4.1. Let uq be a negative rKS potential for which Cp is a positive imaginary number 
for all p — I,... ,P, and let A > 0 be such that the complex critical points ofh{-', d) are all simple 
and correspond to distinct nonzero values o/Re{—//z(z; d)). Then there exists a family of contours 
{Wm]m^i such that the modified system NA^(d)v(d) = Nb^(d) having integrals over contours 
in place o/{C„)^^j is suitable. 


The condition on in Proposition 


4.1 


implies that mq is a negative linear combination 
of Lorentzian profiles 2vp/((x - Upf + vp, Zp — Up + ivp. The proof of Proposition 4.1 can be 
found in Appendix We now show that for suitable systems, the solution v(d) converges to a 
well-defined limit as e —> 0. 


Proposition 4.2. Let A > 0 be given. Suppose that the system ( |14[ ) is equivalent via a matrix N 
to a suitable system NA^(d)v(d) = Nb^(d) having integration contours in its P rows. 

Then 

lim v(d) = v°(d), where (26) 

v“(d) = dlim(z-z„)'p-(z;d) = dRes'p-(z;d), p=l,...,P. (27) 

‘ Z-^Zp Z=Zp 


Proof. We give the proof in the case that all critical points associated with contours {Wm]m=x 
are simple, but the same result is true more generally. Let Xm denote the unique critical point of 
—ih{z'. A) on Wm, which is traversed at the local steepest descent angle 0m, and define a diagonal 
matrix D whose elements are Dmm ■- \h''(xm', V27re. Set A^(d) := DNA^(d) 

and b^ DNb^(d). Therefore 


mp 


bl, = -d 


2716 


^ \ \h"{Xm,A)\ r e 

Jw 


~ih(z;A)lc 


■ dz 


\h"(Xm', d)| ih{x^-X)le -ie, 
2776 


W„ Z- Zp 

e dz. 


f 


(28) 
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so applying the method of steepest descent shows that A^(A) A^^(A) and b^(/l) ^ b^*'(/l) as 

e —» 0, where 

A>°/A):^^— and b^^A)-A. (29) 

Xm - Zp 

The limiting matrix A^°(/l) is of Cauchy type, and it is invertible precisely because the critical 
points associated with the contours are distinct. To show that v*’(/l) given by ( |27l i 

is the (unique) solution of A^*'(/l)v°(/l) = b^°(/l), we may use the Residue Theorem; multiplying 
•HU on the left by A^°(/l) and dividing by -l/d gives 


iy. 




P ResT'-(z;T) 

^ _ 

^ Zp- Xm 


1 

Ini 


'Plz; d) 


c z- 


dz - 1, 


(30) 


where C is a counter-clockwise contour encircling the points (the poles of the integrand) 

and the last equality follows by taking a residue at z = oo where the integrand behaves like 1 /z. 
See im for an alternative approach. □ 

Remark 4.1. When the critical points are simple, A^(/l) — A^°(/l) and b^(/i) — b^®(/l) both have 
complete asymptotic expansions in ascending integer powers of e. It follows that the same is true 
of\{A) — v®(/l). If one or more of the critical points is not simple, fractional powers may appear 
in the expansions. 

Theorem 4.1. Let uq be a rKS potential, and suppose that A > 0 is such that there exists a 
suitable modification of the linear system If A lies outside the bulk (true for all A > 0 if uq 
is a positive rKS potential) then fiA) is exponentially small as e ^ 0. If A lies in the bulk, then 


^ _ ^i lnAix.iA) ^ 

as e ^ 0, where 6+ is defined in ( |22| l and if/+ is defined by 

(■y±(d); d) 


I'T^ (x±(T);T)| 


(31) 


(32) 


with given in ((TS]). The error term in can be written as if the critical points of 

h(z', d) are all simple. 


Proof. By Proposition 


4.2 


the coefficients Vp(A) appearing in ( [T3] l tend to limits v°(d) as e ^ 
0. Thus it remains to analyze the absolutely convergent integral over K of and 

convergent improper integrals over K of for p = 1,... ,P. 

For T 6 \ S, there are no real critical points of /i(z; d), and hence the contour of integration 

in ([T3J may be deformed away from the real axis in the direction of decrease of Re{-i7i(z; d)) 
until the nearest critical point or singular point is reached, which proves the exponential decay of 
y6(d) as e ^ 0 (the distance of the nearest critical/singular point to K determines the exponential 
rate of decay). 

For A e S there are exactly two simple stationary phase points at z = x±(d), so by the method 
of stationary phase, if g(z) - uo(z) or g(z) - (z- Zp)~^, 


1 r 

yllne J-<x 


~ih(z;A)le 


g(z) dz 


gi^/4g-ih-(A)/c 


g(x-(A)) ^ e_ 


-'W4g-i/i+(/l)/e 


g(x+(A)) 


V|/i"(x-(d);d))| 
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V|/i"(x+(d);d)| 


+ 0(e) (33) 















where h^{X) h{x±{A)\ A), and we used the fact that as mq is necessarily a negative rKS potential 

(to have positive d € S) we have sign(MQ(x±(/l))) = ±1. Therefore, 


- ie 


-iMte 





0(e), 


where, using uo(x±(A)) - -A, 


(34) 


p 

Gt(A):^-A-Y, 

p=i 


Vp(4) 

■^±(4) - Zp ’ 


(35) 


The e-dependence in G^(A) enters through the coefficients {vp(/i))^^j; we may replace G^(A) 
witfi Gq(A) (in which Vp(A) is replaced with v^p(A) given by ( |27| )) in ( |34| provided the error term 
is replaced with o(l) (unless the critical points of h are all simple; see Remark [4~T] ). Now using 
we apply the Residue Theorem to obtain 


p 

G^(/l) H-T =-T V Res 

p=\ 


'T-(z;T) 
;c±(/l) - z 


2m Jc x±{A) - z 


(36) 


where C encircles the poles zi,... ,zp positively but excludes the points x±(A). Now 'T“(z; A) 
includes a factor of z - x-(A), so to compute Gg(A) + A the singularity of the integrand at z = 
X-(A) outside of C is removable, so as the integrand behaves like -1/z as z —> oo we obtain 
Gq(A) + a = a or simply Gq(A) - 0. To compute Gg(A) + A we have the same behavior of 
the integrand for large z but now there is a simple pole outside C at z = ■r+(/l), so taking it into 
account gives Gy(/l) + A = A - A'i^^(x+(A); A) or simply Gq(/ 1) = -/l'T“(;ii:+(/l); T). Therefore, 
using may be written as 


m = - 


^ /1|T* (^+(/l); /1)| jMle„-ihXX)lE„-i,l,^U) 

e V|/!"(x+(/l); T)| 


■oie-^'h 


The phase is then simplified using Lemma [TT] and tfie amplitude is simplified by ( [T^ . 


(37) 


□ 


Corollary 4.1 (Matsuno’s modulus formula ifT^ l. Let uq be a negative rKS potential. If A e S 
then 


fime^(/l)r = 27rA(x+(A) - x_(/l)). 
6—»0 


(38) 


Matsuno obtained a result slightly more general than this (he omits the negative rKS condition 
and accordingly generalizes the limit formula) by an argument using trace formulae, identities 
that involve the modulus of fl but not its phase. Theorem 4.1 both makes Matsuno’s result rigor¬ 
ous for a class of potentials uq and also corrects and proves a conjecture in the thesis of Xu 
regarding the phase of p. 


5 . Discrete Spectrum T < 0 

5.7. Evans Function and Eigenvalues 

We study the asymptotic properties of the Evans function ([^ for rKS potentials uq in the limit 
e —» 0 for T < 0; according to (0 it is sufficient to consider only positive rKS potentials. The 
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essence of our strategy is to calculate the determinant in (|^ for e > 0 small by first applying 
the method of steepest descent to the individual entries (|^ of the matrix A(/l), and then compute 
the leading term of the (finite) determinant D(A). However, as in Section]^ it is useful to apply 
this method not to D(A) - det(A(/l)) but rather to D{A) - det(A(/l)) where A(/l) = NA(/l) for 
some invertible N. Otherwise, the computation of the determinant may produce cancellation to 
all orders of the steepest descent expansions of the matrix elements, which only means that the 
dominant contribution to D(A) arises from terms beyond all orders in those expansions i.e., terms 
arising from sub-dominant critical points. Obviously D(A) and D{A) have exactly the same zeros, 
which is all that is important for obtaining the discrete spectrum. Here we choose N so that the 
entries A,„p{A) of A(/l) have exactly the same form as (|^ but with the contours replaced 

with others such that the resulting matrix A(/l) is suitable for asymptotic analysis in 

the sense described in Sectionj^with one exception: the contour of all integrals appearing in at 
most one row of A(/l) will be allowed to pass over either one critical point or two real critical 
points on the zero level o/Re{-//!(z; d)). Such a choice of N has the precise effect of making the 
sub-dominant yet essential contributions to D{A) the dominant contributions to DiA). 

Again, as proof of principle that in a general case a suitable matrix A(/l) = NA(/1) can be 
found, we offer the following. 


Proposition 5.1. Let uq be a positive rKS potential for which Cp is a negative imaginary number 
for all p — I,... ,P, and let A <0 be such that the complex critical points ofh{-', d) are all simple 
and correspond to distinct nonzero values o/Re{-//!(z; d)). Then there exists a family of contours 
such that the modified matrix A(d) NA(d) having integrals over contours {in 
place o/{Cm)^,^[ is suitable and at most one of the contours, say Wq, may pass over either one 
complex critical point or at most two real critical points x±{A) for which Re{—f/!(;r±(d); d)) = 0. 


We prove this proposition in Appendix[B] and therein one can also find an elementary example 
for which N = I, i.e., the original matrix A(d) is already suitable, and an othe r for which N is 
nontrivial. Note that the condition on the coefficients in Proposition 5.1 implies that uq is 


apositive linear combination of Lorentzian profiles 2vpl{{x-UpY+ Vp)w\\e.x& Zp - Up + ivp. We 
now show how to calculate DiA) for a suitable matrix A(d). 


Proposition 5.2. Let uq be a positive rKS potential, and let A < 0 be given. Suppose that for 
some invertible matrix N, the matrix A(d) := NA(d) is suitable. If in addition A lies outside the 
bulk, then 

DiA) = t/,(d)(l + Die)), (39) 

where UYA) is a nonvanishing function defined by while for A in the bulk, 

DiA) = y,(d) (sin(;rr,(d)) + Die)) (40) 


where 14(d) is a nonvanishing function defined by and 

TM) r ix+id) - x-irf)) dp + ^ifi+W - fi-iX)) (41) 

with if±iA) being defined by ( |32| l. The error terms in ( |39| l and (|4gi are uniform on compact 
subintervals of the indicated sets. 

Note that the formula ( [39] l is consistent with the bound (|^ on the discrete spectrum; there can 
be no eigenvalues outside the bulk S, so it is impossible for DiA) (or equivalently 5(d)) to vanish 
there. Note also that at the soft edge d^in = - max{Mo), we have i/'+(dtnm) - lA-(dmin) = 7!’/2. 
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Proof. Let Q be the index of the row of A(/i) whose integration contour may pass through one or 
two real critical points. Since contains exactly one dominant critical point x,„ - x„,(A) 6 
for m ^ Q, the method of steepest descent gives 


Amn(d) — 


2ne 


\h"(Xm, d)! 




■0{e) 


(42) 


as e ^ 0, where hm{A) := h(xm(A); A), and 0,„ is the steepest descent direction with which W,„ 
traverses the critical point x^. This same formula also holds for m - Q, provided d 6 M \ S. 
On the other hand, if d 6 S, the contour Wq traverses two simple real critical points x± - x±(d) 
by the Klaus-Shaw condition, both of which contribute at the leading order to Aqp(A): 


Aqp(A) - 


2ne 


'\h"{x^{A)-A)\ 


-g'W4g-f/!+(/!)/e 


1 


x+( A) - z. 


+ 0(e) 


2ne 


\h"{x-{A)-A)\ 




1 


x-(A) - Zp 


+ 0{e)\, (43) 


where /r±(d) are defined in Lemma 3.1 Aside from the Cauchy factors (x,„ - Zp) ‘ and (x±(d) - 


Zp) * there is no dependence on the column index p in these leading terms. Therefore, by row- 
multilinearity of the determinant, if d 6 \ S, 


5(d) = (2;re) 


pji 




n 


J V|/!"(jCm; d)| 


det(CH-6)(e)), 


(44) 


where C is the P x P Cauchy matrix C^p ■- (^m “ Zp)~^. This matrix is invertible because the 
complex critical points are distinct, as are the poles Hence if the nonvanishing 

function U^{A) is defined by 


U,(A) (.2ne)^/^ 




n 

Vm=l 


y/\h''(xm;A)\^ 


det(C), 


(45) 


the proof of ( [39l l is complete. 

On the other hand, if d e S, then at the cost of a factor (-1)^“^ we may assume that Q - P 
and we then split the last row before computing determinants to obtain 


5(d) = (27re)^/2 


(p-\ 


n 

^m=l 


V|/T'(x„;d)|^ 


^iKlA^-iKWIe jgj (Jet (C') 

+- . -+ 0{e) 


V|/*"(x+(d);d)| 


^|\h"{x_{A)■A)\ 


(46) 


where the matrix C* is the Cauchy matrix C with the last row replaced by ((x±(d) - 
Zi)”',..., (x±(d) - Z/))^5 j and we used the fact that |e~'*±(4)/e| _ j Next we observe that Il24ll 


det(C+) ^^(d) 
det(C-) = -f’(d) 


'I'-(x4d);d) 


x+(d) - X-(d) 
y+(x-(d);d)* 
x+(d) - x_(d) ’ 
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(47) 






























are defined by ([T8]l and where P{A) 0 is given by 


P(A) := 


ni<i<7<r'-l(-^; -^i) ni<i<7</>fc Zy) 

Y\m=l Yln^li^m - Zp) 


Using ([T^ and Lemma 3.1 noting that because uq is positive Klaus-Shaw 


X oo ^0 

t/oW dx - I (x+(r/) - X-(r/)) dr} 
oo Ja 


-f: 


max|Mo} 


{x+{p) - X-iri)) dr], 


and defining (restoring the factor (-1)^ 


(48) 


(49) 


V.(d) 


-2ii-lf-Q{27:ef/^P{A) 

V-d(x+(/l) - .x_(d)) 


fp-i 


n 

^W=l 






then establishes (|40li and completes the proof. 


(50) 

□ 


Remark 5.1. It is easy to see that if the fixed value A e S is replaced with an e-dependent value 
df with the property that for some fixed A e S and some constant C > 0, \A — dfl < Ce holds 
for all e > 0 sufficiently small, then the asymptotic behavior of the matrix elements A„,p(A^) can 
also be calculated by the method of steepest descent, with results similar to One need 

only replace the functions h,„{A) — hix^iAf A) and h±(A) — h(x±(A);A) in the exponents on the 
right-hand sides of (|4gi-(|43 by h{Xm{A)', /If) and h(x±(A); /If) respectively. 

Proposition |5.2| allows us to establish a number of important asymptotic properties of the 
discrete eigenvalues for positive rKS potentials. 


Theorem 5.1 (Uniform approximation of eigenvalues). Let uq be a positive rKS potential and 
suppose that for each A e K <z S, K compact, a suitable modification A(/l) of A(A) can be found. 
Then there is a constant Ck such that for each eigenvalue A in K there exists /lo < 0 satisfying 
7’f(/lo) e Z, such that |/l — /lol < Ckc^ holds for all sufficiently small e > 0. Likewise, for each 
/!() 6 K satisfying T^Aq) 6 Z there is an eigenvalue A such that the same estimate holds true. 


5.2 


in the 


Proof. The eigenvalues are characterized exactly by D{A) - 0, or using Proposition 
case A e S, sin(7r7’f(/l)) = 0(e), the error term being uniform for A e K c S. Solving this 
equation for T^A) and multiplying by e gives, for some n e Z, 


1 

271 


r 


(x+(rj) - X-(T]))dt] + —e(fi+(A) - fj-i^)) = en + 0(e^) 
2n 


(51) 


uniformly on K. Since the left-hand side is differentiable with respect to A with a derivative that 
is strictly positive on K, the result follows from the Implicit Function Theorem. □ 


Corollary 5.1 (Local approximation of eigenvalues). Fix a closed interval K c S and a positive 
integer J. Under the hypotheses of Theorem \5. 1\ for each A € K the 27-1-1 eigenvalues A closest 
to A are given by 


/I = A -H 2ne 


j + [rf(A)] - rf(A) \ 

x+(A) - X-(A) I 


0{e\ ; = -7,-(7-l),...,7-l,7, 


(52) 
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in the limit e ^ 0, where the error term depends on J and K, and where [■] denotes the nearest 
integer function. Note that |[7’f(A)] — 7’f(A)| < 1/2. 

Proof. We set /i = A + A and write the exact eigenvalue condition as eT^{A + A) - en + 0{e^) 
for some integer n. Writing the integer n in the form n - [^^(A)] + j for an integer j e [-7, J] 
and applying Taylor expansion, 

e7’'(A)A + 6>(A2) = e (/ + [^.(A)] - ^.(A)) + Oie\ (53) 


because derivatives of are bounded uniformly for A e K. Now we seek solutions A = d - A 
that are 0(e), so the error terms may be combined and we may solve for A; 


A-A 


e(j + [r,(A)] - r,(A)) + 0(e^) 

er'(A) 


(54) 


Since eT'^(A) - (v+(A) - x-(N))!(In) + 0(e) uniformly for A e K the proof is complete. □ 


This result shows that the eigenvalues are locally equally spaced with 0(e) spacing that de¬ 
pends on the point A of local expansion. Ignoring the details of the equal spacing and the offset 
of the grid given by the term [^^(A)] - T^(A), we reproduce a result obtained by Matsuno m 
by formal asymptotic analysis of trace formulae (conservation laws). 

Corollary 5.2 (Matsuno’s density formula ifTSll l. The asymptotic density of eigenvalues at a 
point A in the bulk is pyi(A)le, where Matsuno’s density is 

Pm(X) := ;^(.s:+(d) - V-(d)). (55) 

ZTT 

In other words, the number N[a, b] of eigenvalues in a subinterval [a, b] of the bulk tB satisfies 

1 r* 

N[a,b]-- I pM('^)dA +0(1), e ^ 0. (56) 

e Ja 

For all positive rKS potentials uq, Matsuno’s density pM vanishes at the soft edge and blows 
up at the hard edge d = 0. This should be contrasted with the better-known Weyl density of 
eigenvalues for semiclassical 1-D Schrodinger operators, which is typically hnite and nonzero at 
both hard and soft edges. 

Remark 5.2. In 1 1771/ it was shown that if uq(x) — 2/(1 -i- x^), then for e — I IN for a positive 
integer N, the eigenvalues are precisely the roots of the equation Lf^l—INA) — 0 where Lfj is the 
Laguerre polynomial of degree N. In this case, the distribution of eigenvalues is well-known in 
the theory of orthogonal polynomials and random matrix theory. In the latter context Matsuno’s 
density for the case uq(x) — 2/(1 + x^) is known as the Marcenko-Pastur law am ; from this 
context that we adapt the terminology of soft and hard edges of the (discrete) spectrum S. From 
the formula 01 one can see that for this potential, when e — l/N the reflection coefficient f 
vanishes identically. 
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5.2. Phase Constants 

Fix a number A e S and for each e > 0 sufficiently small let be the eigenvalue closest to A, 
choosing df > A if there are two equidistant. According to Corollary [AT] there is some constant 
C > 0 depending only on A, such that |A- dfl < Ce holds for all e > 0 sufficiently small. Denote 
by and Of(v) the phase constant and normalized eigenfunction, respectively, associated to the 
eigenvalue A^. Then (|T5|) takes the form 




InAf 


= lim 

R^+oo 


- lim 

i?—»+oo 


X oo 

0^(x)* (xcl)^(v) -l)dx 

CO 

X R ^ nR 

R 


£ 


xIJx) dx - 

2nA, 

ie 

xh{x) dx - — 

R ^^6 


J 

J-R 


<t>,(;r)* dx 


(57) 


h(x) 


2tiA, 




In the last step we deformed the contour [-R,R] to a semicircle in the lower half-plane where 
0(z*)* is analytic and satisfies <l>(z*)* = + 0(z^^) as z ^ oo. 


Proposition 5.3. Let A e S be such that the complex critical points ofh(-;A) are simple, and 
set L — L(A) 1 -H max{|v+(A)|, |v_(A)|). The function fix) defined for x E in •nzl has the 
following properties. 

1. There is a constant Kq — Kq(_A) such that 


fix) - 


e 

2nA^x^ 


< ^ 
|jcp ’ 


\x\ > L 


(58) 


holds for all e > 0 sufficiently small. 

2. There exists a constant Ky - KfA) such that 


\hix)\<Ku \x\<L 


(59) 


holds for all e > 0 sufficiently small. 

3. For each real x + x±(K), 


lim/f(.r) = /o(x) 
6-^0 


/f'(x_(A),j;+(A))(-^) 

V+(A) - V_(A) ’ 


where X{a,b){x) denotes the characteristic function of (a,b). 


(60) 


Remark 5.3. Comparing with the definition of fix) given in shows that for v e M 
the modulus of the normalized eigenfunction d>f(x) behaves, in the limit e ^ 0, like a large (of 
size multiple of the characteristic function of the interval bounded by the turning points 

v±(A). This asymptotic behavior can be compared with the more familiar result, based on the 
WKB method, that eigenfunctions of a Sturm-Liouville problem, while also being asymptotically 
confined to a “classically allowed” interval (a_(A), a+(A)), nonetheless have a decidedly non¬ 
constant amplitude on this interval, diverging at the endpoints typically like \x — a:±(A)|“^^'*. This 
draws a strong contrast between the nonlocal BO scattering problem and scattering problems 
for many other integrable equations like KdV involving spectral analysis of purely differential 
operators. 
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Proposition |5.3| is proved in Appendix [C] Its main purpose is to prove the following theorem. 
Theorem 5.2. For each A e S, 


limy^ = -x(A), x(A) ^ (x+(A) + X-(A)), 

6-^0 Z 


( 61 ) 


where is the phase constant for the eigenvalue nearest A. 

Proof. Since df ^ A < 0 as e ^ 0, from we see that Im{yf) ^ 0 as e ^ 0, so it remains to 
analyze the real part. Next, recalling L - L(A) as defined in Proposition |5.3| and using property 1 
from the same yields, for each R > L, 


dx 


r 


r edx 


6 

( xl^(x) dx 

M 

< 

jjR 27rAfX 

_ 

+ 



L 

X 


|x| f(x) ■ 
Kgedx 


2nAtX^ 


dx 


(62) 


'tR X^ 
L 

2Koe 


:= [-R,-L] U [L,R], 


an upper bound that is independent of R > L and that tends to zero with e. Hence 


lim 7 ^ = lim I 

e —»0 6-^0 


xf{x) dx. 


(63) 


so applying the Lebesgue Dominated Convergence Theorem to calculate the latter limit, making 
use of properties 2 (integrable e-independent domination) and 3 (pointwise limit), gives 


limy, 


nt- nJf+A) 

- I xlo(x) dx - - I 
J-L d.t_(A) 


xdx 


xXA) X+(A)-X_(A) 


Evaluating the integral then proves 


(64) 

□ 


Remark 5.4. Theorem 5.2 shows that the limit ofy^ is a purely real number. On the other hand, 
Imjyf) = —iel(2Af) cannot always be neglected outright. For example, it plays a pivotal role in 
the small-dispersion analysis of certain determinantal T-functions arising in inverse-scattering 
theory l\20\l . 

A version of the asymptotic formula ( |M] l characterizing the real part of y^ for small e was 
hypothesized, based on asymptotic analysis of the BO inverse scattering problem, in Il20l . The¬ 
orem [5|^provides a careful statement and is the first rigorous result on the asymptotic behavior 
of the phase constants, which (like the phase of the reflection coefficient) cannot be captured via 
trace formulae. 


6. Numerical Verification 

We illustrate the accuracy of our asymptotic formulae for the scattering data by comparing 
them to exact calculations in the case of the rKS potential: 


uo(x) - -vi 


1 


x-i X - (/ H- 1) 
16 


-I- c.c. 


(65) 



















for V = +1. The graph of mq is plotted in Figure]^ for v = 1, which confirms the Klaus-Shaw 
condition. This rKS potential is positive when v = 1 (negative when v = -1). The bulk S 
consists of the interval (-max{Mo),0) when v = 1 with max{Mo) ~ 5.07308 and (0,-min{Mo)) 
when V = -1 with min{Mo) ~ -5.07308; see Section]^ and Definition |3.2| We selected this 
rKS potential so that it is not even about any point, a property making the phase constants jj 
nontrivial to calculate and interesting to compare with small-dispersion asymptotics. 

We compute the exact eigenvalues Aj for the potential ( [hS] ) with v = 1 using the Evans function 
(H). To simplify computations only values of e > 0 for which icije 6 Z and /ca/e 6 Z are 
considered. Since ic\ and ic 2 are integers, this requires e = 1/m for some m e N. In this case, the 
integrals Q defining the elements of the matrix A(/l) can be calculated explicitly by the Residue 
Theorem and the exact eigenvalues are thus obtained as the roots of a polynomial; see HD. 

In Figure]^ we show the exact eigenvalues (black dots) and compare them (i) with their uni¬ 
form approximations obtained from Theorem |5.1| by solving T^{A) - n for positive integers n 
(overlaid squares) and (ii) with their local equally-spaced approximations described by ( |5^ in 
Corollary |5. 1 1 (overlaid circles), under the scaling (d - A)/e. As expected, the uniform approxi¬ 
mation’s squares track the black dots very well, while the circles do so best near the point d = A. 
In Figure 1^ we present histograms of the exact eigenvalues to highlight their distribution on M . 
The histograms clearly match better and better, in the limit e —> 0, the density pM(d) of Matsuno 
(Corollary|5.2|i, here normalized to have integral equal to unity. 


e e 



Figure 3: Exact eigenvalues for the rKS potential with V = 1 as a function of e (Black dots) with the uniform 
approximations (squares) and the local approximations based at various values of A 6 S (circles) overlaid. Note that the 
horizontal axis is the rescaled local coordinate (A - A)/6 with A = -4 (top left), A = -3 (top right), A = -2 (bottom 
left), and A = -1 (bottom right). 


In Figure]^ we illustrate Theorem |5.2| by plotting the real parts of the three exact phase con¬ 
stants corresponding to eigenvalues d^ closest to A = -1,-4, -5 in the bulk S. The limiting 
values predicted by Theorem |5.2| are indicated with dashed lines. The exact values of Re{ye) 
were computed using the alternate formula for y^ presented in lfT9ll . 

The reflection coefficient is calculated exactly from ( [T3] l with the coefficients Vp(d) computed 
by solving the linear system Q for v = +1. To aid in numerical computation we deform the path 
of integration in ( [T3] l into the complex plane to exploit the exponential decay of the integrand. 
We compare the exact reflection coefficient [i with its asymptotic approximation (Theorem |4.1| i in 
Figures [^and[^ In Figurej^we plot the normalized magnitude ^[e\|3\ as a function of the spectral 
parameter d. The left panel (v = -1) shows that as e —> 0, yfe\l3\ indeed approaches the expected 

17 














1.5 

e = 1/16 

1.5 

e = 1/32 

1 

48 eigenvalues 

1 

96 eigenvalues 

18 bins 

27 bins 

0.5 

0 


0.5 

0 



1.5 

1 

0.5 

0 


A 


A 


e = 1/64 

1.5 

e = 1/128 

192 eigenvalues 

1 

384 eigenvalues 

39 bins 

56 bins 


0.5 

0 



Figure 4: Histograms of exact eigenvalues for the rKS potential with V = 1 illustrating the distribution of eigenvalues. 
The limiting curve (solid black) is Matsuno s density pm(^) normalized to unit mass. 



Figure 5: Three values of Rely^} (solid lines) as a function of e for the rKS potential j65) with v = 1. The limiting 
values (dashed lines) are from |6T) evaluated for A = -l,-4,-5. The values Rely^} are computed at each 6 with the 
exact eigenvalue chosen to be the closest to A. 


limit — whose support is the bulk — given by Matsuno’s modulus formula (Corollary |4.1| l. The 
right panel (v = 1) shows that ■\/e\/5\ —> 0 for positive rKS potentials as predicted by Theorem |4.l| 
In Figure |7]we plot the derivative of the phase of yS (computed indirectly using elm{j8'(/l)/j0(/l))) 
and compare with the corresponding limiting curve predicted by Theorem 4.1 namely - - 

-x+(A). 
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Appendices 

A. Proof of Proposition |4.1 1 

Since Im{cp) > 0 for each p, the integrands for the elements of the matrix A(/l) and the vector 
b(/l) defined for T < 0 by are integrable at all of the singular points {Zp}^^^ in C^, so when 
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Figure 6: Left: V^|/?| as a function of A for the rKS potential {65) with y = -1 for e = 2, 1, 1/2, 1/4, 1/8, 1/16 (solid 
black curves). The apparent limiting curve (dashed-blue) is obtained from Corollary |4.1| For 6 and A both small the 
graphs become difficult to compute and are not plotted. Right: Same as the left panel but for v = 1 and 6 = 4, 2, 7/4, 
13/8, 3/2, showing convergence to zero. 





Figure 7: Plots of the rescaled derivative of the phase for the function /? as a function A. The limiting curve -x+(/l) 
(dashed-blue) is obtained from taking the derivative of the leading order phase -Q+iA) in ID- 


/i < 0, in each row for which the contour is C,„ = (the general case) we may write the 
integrals in Q in the form 

r f (A.l) 

Jc,„ 

where either g(z) - (z-Zp)"^ or g(z) - 1; see Figurej^ Note that the term with n - m in the sum in 
necessarily has a nonzero coefficient because otherwise Cm = £o(Zm) as icmie is a negative 
integer (the special case). Carrying out the analytic continuation of A(/l) and b(/l) through 
and taking the boundary value on we arrive at formulae for the matrix elements of A^(/l) and 
b^(/l) for d > 0, in which the infinite tails of all of the contours €„-\{Zn) for n - 1,..., m or toizm) 
are all rotated through the left-half z-plane so that they now originate at -ioa to the left of the line 
Re{z) = Re{zi), and at the same time the integrand is analytically continued as necessary. We 
define an intermediate set of contours such that Fi connects -ioa (left of Re{z) = Re{zi)) 

to zi and then for m - 2,... ,P, Ym connects Zm-i to Zm, with in the domain of f(z). Then the 
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rotation of the integral over (n-xi^n) in ( |A.1[ ) may be written as a linear combination of integrals 
with the same integrand over contours Fi,..., with the coefficient of the integral over y„ 
being 1. The same holds for the integral over £o(Zm) that arises in the special case. The other 
coefficients in these combinations are nonzero factors arising from analytic continuation and the 
multi-valuedness of f(z). It is then clear that there is a lower triangular matrix M with nonzero 
diagonal entries M,„m = (1 - in the general case and M,„„, = 1 in the special case such 

that for m - 1,..., f*, for T > 0 we have 



(A.2) 


Next we show how to exchange the contours for a different set more 

amenable to asymptotic analysis and defined by the following procedure. Consider the level 
curves of R(z) := Re{-ih{z', T)) = L for levels L increasing from L - -oo to a sufficiently positive 
value. Under the assumption that Re{cp) = 0 for all p, R(z) is a harmonic function in the domain 
C \ {zi, ■ ■ ■ ,zp,z’l, ■ ■ ■ ,z*p]. The landscape of R{z) can be visualized as being based upon a plane 
with slope T > 0 in the positive imaginary direction; this is the contribution of the term Az in 
h{z', T). Superimposed on the plane are P infinite mountain peaks at points z*,..., zj, in the lower 
half-plane and an abyss consisting of bottomless pits at the corresponding points z\,. ■ ■ ,Zp in 
the upper half-plane. The latter are all features contributed from the term /(z) in h{z', d). If we 
visualize the level L as the sea-level on the landscape of R{z), then when L is large and negative 
there will be small distinct lakes formed in each of the P pits and otherwise the upper half-plane 
will be dry land. Under a steady rain the sea level L will rise, and at the critical values of L = Riz) 
corresponding to the critical points of h with positive imaginary parts, some lakes will fuse; the 
number of lakes thus decreases in steps as L increases. Meanwhile, the tide steadily rises from 
the lower half-plane, surrounding the mountains and forming islands that become isolated from 
each other at critical values corresponding to the critical points of h with negative imaginary 
parts. Eventually the lakes in the upper half-plane will begin to merge with the ocean as well 
as with each other, and ultimately (for large enough L) the shoreline in the upper half-plane will 
consist of a single curve stretching from horizon (Re{z) = -oo) to horizon (Re{z) = -i-oo) as all 
lakes will have been subsumed by the rising tide. Note that when L - 0 the entire real axis is a 
component of the level curve. 

Since the values of R(z) at its critical points in are distinct and nonzero, as the sea-level 
L increases through a critical value, exactly two lakes will fuse, or a single lake will fuse with 
the ocean, with the point of fusion in each case being a critical point of R. If L = 0 is a critical 
value of R, then there are at most two real critical points x±{A), distinct for A e S, and at L = 0 
these are distinct fusion points where a single lake meets the ocean. At the lowest critical level 
L, the fusing lakes each contain exactly one of the pits (or none, in the special case of 

fusion with the ocean), and we define a contour W connecting the two pits (or connecting one 
pit with the point -/oo in the ocean) and lying below sea-level except at the critical point which 
it crosses locally in the steepest descent direction (if the lowest critical level is L = 0 we make 
the connection via X-(A) 6 M). Increasing L then completely encloses W within the new lake or 
ocean formed by the fusion. Suppose that at some non-critical level L, each lake contains a chain 
of contours W each of which connects a pair of pits (or for the ocean, exactly one of the links 
W connects a pit with -ioo) such that all enclosed pits are vertices of the chain. This inductive 
hypothesis is certainly true for levels L just above the lowest critical level L. Increasing L to 
the next higher critical level, define a new contour W to pass over the fusion point locally in 
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the steepest descent direction (choosing x_(/i) as the fusion point if L = 0) but otherwise lying 
below the sea-level and connecting one of the finite endpoints of the chain (a pit) in each fusing 
lake or ocean. Thus two chains are joined into a longer chain that is enclosed within the formed 
lake/ocean for slightly larger L, and the induction argument closes. 

Once all of the lakes have been subsumed by the ocean, we have assembled a single chain 
of contours W with one finite endpoint (a pit), one infinite endpoint at -/oo, and whose internal 
vertices are all of the remaining pits. We now label the links of the chain in order from the 
infinite end to the finite end as VTi,..., Wp, and we orient each link consistent with this labeling. 
Each contour W,„ contains exactly one critical point z - x,,, of h{z', d), the global maximizer of 
R restricted to W^, with x„, 6 unless both m - 1 and A e S, in which case xi - X-(A) e M. 
We concretely define the multi-valued function h{z', d) on W„, by arbitrarily selecting a branch at 
z = x,„ and analytically continuing along W^. The last step is to observe that we may express 
integrals of along the contours as linear combinations of integrals of the 

same integrand, with the branch chosen as indicated above, along the contours and 

vice-versa. But this is easy, because both chains of contours and span the same 

set of nodes {Zp}p^Q, Zo := -loo, so taking into account possible monodromy of the integrand, 
there indeed exists an invertible matrix K such that 

f y r (A.3) 

Aw„ ^ Jy„ 

Combining our results, we see that with N KM“' the modified system NA^(d)v(d) = 

Nb^(d) having contours in place of is indeed suitable as claimed. 

B. Proof of Proposition 

We prove Proposition |5.1| in the most interesting case that A e S and hence there exist two real 
critical points x±(d), first associating a certain tree graph to the level sets of R{z) Re{-ih(z', d)) 
generated by the potential uq. The description of the level sets R{z) = L is analogous to that 
presented in Appendix but the landscape is inverted, sloping downward as Im{z) increases, 
and with “pits” and “mountains” interchanged. 

B.l. Construction of a Tree 

To construct a unique tree associated to the initial condition ([^, we consider the level sets 
R{z) = L as L is varied. For large L > 0, the level sets of R in will consist of P disjoint loops 
each bounding a single “island” containing exactly one of the poles each of which is an 

infinite mountain peak since icpfe > 0, and each of which is associated with a distinct leaf of the 
tree. The P islands will be associated with edges in the tree terminating at the leaves. As L is 
decreased to the first (largest) critical value, two islands will fuse at a critical point of the function 
Riz) to form a single island containing two of the points as L is decreases slightly further. 

We associate these fusion points of islands with tree nodes having two outward-directed edges 
representing the two fusing islands and one inward-directed edge representing the newly-formed 
larger island. We include an additional node with the label x± for the island that merges onto the 
real lin^(part of the level set L = 0); for L < 0 a “continent” invades that may fuse with 


’ Since the potential is assumed to be Klaus-Shaw, the level sets (islands) merge onto the real line only at the points 
x±. 
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islands if there are any negative critical values of R. For convenience, we include an additional 
edge in the tree connecting the node associated with the lowest critical value L to a root node 
labeled “-oa” at height L - -oo. Finally, we label the leaves of the tree . ,(p in order of 
decreasing “height” L of their parent nodes (breaking ties arbitrarily) and we label each node of 
the tree Xp to correspond with its associated critical point. The labeling of the leaves induces 
a permutation matrix P mapping the indices of the poles to those of the corresponding 

leaves. We have thus constructed a rooted tree with P internal nodes, P leaves at height 

L - oo, and a root node -oo at height L - -oo; see tree A in Figure [BT^ for an example. 

B.2. Tree Pruning Algorithm 

We obtain suitable contours from the tree graph by using a “pruning” algorithm to 

identify all but one of the leaves of the tree (associated to poles via P) with corresponding 

nodes (associated to the critical points In a manner to be explained, it is thus possible to 

associate a critical point and hence a level set which will be partly enclosed by the contour W„,; 
see the bottom panels of Figures [BT9| and [B . 1 0| 

The first step is to chop (divide) the tree at the node x±, with the split node x± becoming 
simultaneously a leaf of one subtree and the root of the other. We are thus left with two trees: 
one having Q leaves, Q - \ internal nodes, and a root node labeled x±, and another having 
P- Q+\ leaves (one of which is labeled x±), P - Q internal nodes, and a root node labeled -oo. 
The leaves of the subtree with root x± are necessarily those labeled . ,{q (all descendants of 
an internal node from a level L > 0) while the leaves of the other subtree are x± and {q+i, .. .,ip. 
The tree-pruning algorithm for the subtree with root x± is the following. 

For m - 1,..., Q - 1 

• Prune the leaf and remove its incident edge and parent node (joining the two remaining 
edges incident on the removed node). 

• Define as the formal sum of leaves 4 already pruned from the subtree that were originally 
descendants of the removed node and associate to the removed node (and hence to a 
critical point x,„). 

Note that for each m = . ,Q - \ the sum contains and possibly some of the leaves 

fi,..., fm-i, but no leaves with indices exceeding m, and that only the leaf {q does not appear in 
any sum. We now repeat the algorithm on the subtree with root -oo, with the index m running 
over the range Q+ 1,..., P. Again the sums i„, will have a lower-triangular structure with 
appearing in the sum Sm, but the leaf x± will not appear in any sum. After the algorithm is applied 
to both subtrees, we set sq = H- ■ ■ ■ H- {p and associate it with the node x±. 

In this manner, the algorithm associates the formal sums uniquely to the complex 

critical points {jCplpZj. Moreover, the coefficients in the formal sums form a P x P zero- 

one matrix K with block structure: 


Lq-i 

0(e-i)x(p-e+i) 

O(/>-0xe 

Lp-e 


(B.l) 


where L,q-i and L,p-Q are square lower-triangular matrices of the indicated dimension with ones 
on the diagonal, and where the dots stand for a single row of P ones: 1^. The matrix K is 
obviously invertible. 
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We now relate the formal sums to contours in the complex plane as follows. First 

consider m Q, and let x,„ be the complex critical point associated with s^', the level curve Riz) - 
R{xm) = L is therefore a “figure-eight” shoreline of two fusing islands (or one island fusing with 
the continent), and one lobe of the figure-eight contains exactly the poles Zp corresponding to 
the leaves in the formal sum Sm- We then take Wm to be the contour consisting of (i) a short 
arc of the steepest descent path for R(z) over the critical point Xm and descending to the level 
L - 5 for some small d > 0, (ii) arcs of the level curve R{z) - L - 6 meeting the arc from (i) 
and ending on one or two branch cuts of / emanating from singularities associated to leaves in 
the formal sum i„„ and (iii) two vertical rays along those branch cuts going to /oo. In the case 
that the figure-eight has an unbounded lobe (the continent), the algorithm guarantees that the 
bounded lobe will always be selected. So, the contour Wm thus contains a unique critical point 
of h, which is the maximizer of R on the contour. The exceptional contour Wq is then related 
to the corresponding formal sum sq - . ,{phy being defined to be homotopic to Up (thus 

enclosing all of the singularities {Zp)^^j) and to pass over the real critical points x±{A) such that 
they act as simultaneous maximizers. 

Finally, let Up, p - 1,..., P, be a U-shaped contour like those in Figure [^but enclosing only 
the point Zp along with its branch cut. The obvious relation C,„ = Um - U\ + ■ ■ ■ + U,n for m - 
I,... ,P clearly defines a lower-triangular invertible zero-one matrix M. The contours {Wm]m=i 
therefore appear in the rows of the matrix A(/l) = NA(/l), where N ;= KPM^' is invertible, and 
as each contour Wm contains exactly one dominating critical point, with the exception of Wq 
which contains two both at level L = 0, the proof is complete. 


A B C D E 

L ^2 ^3 ^4 ^5 ^2 ^3 ^4 ^5 -^3 ^4 ^5 ^4 ^5 ■^4 



Figure B.8: Example tree (A) and the tree pruning algorithm (B-E). Note that tree A has been chopped at the node 
x± which is “shared” by the subtrees. The pruning algorithm on the subtree with root x± gives rise to the following 
assignments: X 2 i-» (^i = ([), X 4 i-* (i 2 = (i + £ 2 ), and xi i-* (^3 = £ 3 ). Similarly, the pruning algorithm on the subtree 
with root -00 gives rise to the assignment X 3 i-* (is = £ 3 ). Finally, i-» (ia = £[ + £2 + £3 + £4 + £ 3 ). 


B.3. Elementary Examples 

Consider the positive rKS potential mq with P - 3 and data 

(zi,Z2,Z3) - 5 (21,2 + 41,-2 + i) 

, (B.2) 

(Ci,C2,C3) = - 51 ( 2 , 1,4), 

and let T = -1. While it is not obvious from the data ( |B.2[ i, from a graph of uq one can see that it 
satisfies the Klaus-Shaw condition. The scheme for selection of appropriate contours [Wmil^^y in 
this case is illustrated in Figure [B^ Obviously Wj is homotopic on the domain of analyticity of 

23 ' 










/ to Cj = UJ (see Figure j^, so here N = I and hence A(/l) = A(/i) and D(A) - D(A). Note that 
as desired, only Wp =3 traverses more than one critical point at the same level of Re{-ih(z', -1)). 




Figure B.9: Top left panel: for the P = 3 positive rKS potential uq with data given by jB.2) , branch cuts of / emanating 
from the branch points {zi ,Z 2 ,Z 3 i are shown with zigzagged lines, the solid curves are the levels of Re{-ih(x', -1)}, and 
the dashed curves are the levels of Im{-// 2 (A:;-1)} (steepest descent/ascent). The intersection points mark the critical 
points of h(x;-l) which are numbered according to the values RQ{-ih{xj',-\)) except for those {x±) on the real axis. 
The domain Re{-ih{z', -1)] < 0 is shaded. Top right panel: the tree associated with A(jc;-1), with the permutation P 
shown in the inset. Bottom panels: the contours = t/^ = Cm for which the integrals Ampi-l) are exponentially 
dominated by a contribution from neighborhoods of the critical point(s) over which the contour passes. The level curve 
Re{-?/i(z;-1)} = const, containing the traversed critical point(s) is plotted, and the domain Re{-z/i(z;-1)} < const, is 
shaded in each case. 


To see that such a trivial outcome is not generally the case, we consider now replacing 
with 


{zuZ2,Z3) - 5 (-4 + 6i, 2i, 2 + i) 
(ci,C2,C3) = 1,3), 


again yielding a positive rKS potential with P = 3. See Figure B.IO 
the algorithm is 


The matrix N returned by 



1 

0 

0 

0 

0 

1 

1 

0 

0 

-1 

0 

-1 

1 

KPM ' = 

1 

1 

1 

0 

1 

0 

1 

1 

0 

= 

0 

0 

1 


0 

0 

1 

1 

0 

0 

1 

1 

1 


1 

0 

0 


(B.4) 


C. Proof of Proposition |5]^ 

Recall the vector <f>{A) characterized for an eigenvalue 4 < 0 by the rank P - 1 homogeneous 
system A(/1)0(/1) = 0 and the normalization condition \^^{A) - A, where 1 ^ := (1,1,..., 1). 
Recall also the eigenvalue A^ nearest A e S. 


Lemma C.l. For each A e S, the elements cf>p(A^), p — I,... ,P of the vector 4>(A() satisfy 


lim fpiAf) 

6^0 


.0 'P-(z;A) 

Pp := ARes- 


z=Zp Z - X_(A) 
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nfc' ( 

Zp - xfA)) 

Wi*p 

{zp - Zl] 

1 


(C.l) 

































































Figure B.IO: Same as Figure[B^but for the positive rKS potential with jB.3) . 


where is defined by 

Proof. By taking appropriate linear combinations of the rows and multiplying on the left by a 
suitable invertible diagonal matrix D(/Ie) to be determined, we replace the system A(A^)^(Ac) - 
0 by the system D(/lf)A(/lf)0(/lf) = 0, maintaining the rank as P - 1. Next, we replace the 
(redundant) P* row of D(/lf)A(/lf) with 1^, resulting in a matrix we denote as A(/lf). Since 
e(fi(Air) = Ag, (fi(Ag) is a solution of the square inhomogeneous linear system A(/lE)^(/lf) = 
b(/l(r) (0,0,..., 0, /If)^. We prove the lemma by showing that for small enough e this system 

has full rank P and by explicitly constructing the unique limiting solution. 

Recalling ( |42) and Remark [5T| we see that upon choosing 

DM.) := ^ 0 (C.2) 

V 2716 


for m = 1,..., P - 1, we obtain limf^o A(/lf) = A°, where 


- 

^mp - 


.^m(A) 


OT=1,...,P-1, p — I, . . . , P 


(C.3) 


while the last row of A° is 1^ (unchanged in the limit). Also limf^ob)/!^) = b° 

(0,0,..., 0, A)^. By a direct computation one sees that det(A°) + 0 because the critical points 
.ri(A),... ,xp-i(A) are simple and hence distinct, as are the poles z\,.. ■,Zp. 

It therefore only remains to show that the unique solution of the limiting system A°0° = b° 
is given by O’ a fact that is a simple consequence of the Residue Theorem; see Il24l for a 
complementary approach. Indeed, for7«=1,...,P-1, 



p=\ 




1 „ 'B-fcA) 

-Res- 

Zp - Xm{A) z=Zp Z - X_(A) 


A 

2ni 


i 


'P-(z;A) 

(z-xUA))(z-xM)) 
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where C is a counter-clockwise contour about the poles of the integrand On the other 

hand, this integral vanishes because the integrand is analytic in the exterior of C and is 0{z~'^) as 
oa. Similarly, the last row of the product is 




A r '¥-(z;A) 
Jc z- x-(A) 


dz - A, 


(C.5) 


the integral being computed by a residue at z - oo. □ 

Recall that <l>f(;c) = <[)(x; df) is given by with A - df. Also, since d^ is an eigenvalue, 
Of(x) —» 0 as X —» H-oo, which implies that the lower limit of integration in ( [TOl i can be replaced 
by -Hoo when d = d^. Using the definition of /f(x) given in ( |57l i then shows that (the two signs + 
are arbitrary but equal) 


p p 


(x) = Y '^'^4>p(A^)^g(A^yapJx)al^(xT, 


p^l q=l 


where 




1 r 


I 


z-Zn 


h^(z) h(z; dj. 


We prove properties 1-3 of Proposition |5 .3 1 in turn. 


(C.6) 


(C.7) 


Property 1. To prove ( |58| ), suppose that x < -L. Integrating by parts using (d^ -i- uo(z))e 
gives 


a;„(x)^i^- 


-ih^{x)je 


(df + uq{x)){x - Zp) 


e-‘‘'‘^^^l‘u'^{z)dz 
(df H- uq(z))Hz - Zp) 
To control the last term we integrate by parts again: 


/' 


f 


e-ihAz)/edz 


(df -I- uq(z))(z - Zp)^ 


(C.8) 


x: 


e 


~ihe{z)le 


dz 


_oo (df + uo(z))(z - Zp)^ 


- le 


-ihe{x)je 


X 


(df -H uo(x))Hx - Zp)^ 

‘ le-'’^‘^^l^u'^{z)dz 

„ {A^ + UQ{z)f{z- Zp)^ 


f 


2e-i>'Az)le Y 


(df H- Uo{z))Hz - Zpf 


(C.9) 


The integrals in ( |C.8| )-( |C9| l are, unlike that on the right-hand side of ( |C.7| i, absolutely conver¬ 
gent. Since d^ -H uo(x) is bounded away from zero uniformly for e sufficiently small and x < -L, 
and since there exists K > 0 such that |mUx)| < K\xy^ holds for all x 6 M, the dominant term in 


Qp ^(x) comes from the first term on the right-hand side of ( |C.8| l, and we easily obtain 


- i 


~ihe(x)le 


2ti dfX 


A'Vi 


X < -L, p - 1,..., R. 


(C.IO) 


for some constant K' > 0. Using ( |C.10| l in ( |C.6| l and taking into account 1^0(df) = Ac as well as 
the boundedness of ^(d^) in the limit e ^ 0 implied by Lemma C.l proves ( [581 ) for x < -L. For 
x> L, one simply uses (|C.6|l with ai 


-p c(x) in place of Op ^x). 
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Property 2. To prove ( |59l ), it follows from ( |C.6| l and Lemma C. 1 that it is sufficient to uniformly 
bound for < T(A) and ^(x) for x > x(A), because ^ A < 0 as e ^ 0. We will 

prove that \a^ ^(x)| < K holds on ;<: < x(A) for all p - 1,... ,P and all e > 0 sufficiently small, 
with the complementary estimate of |a^^(x)| on x > x(A) being similar. Moreover, from ( |C.10| l 
we see that |a^^(x)| < K certainly holds for x < -L, so it is enough to consider the difference 


Op ^ix) - Op ^i-L) for -L < X < x(A). 


The limiting exponent function h{z'. A) has a unique critical point at z = x_(A) in the interval 
(-L, x(A)) at which /!"(x_(A); A) = Mg(x_(A)) > 0 because uq is positive Klaus-Shaw. The 
equation h{z'. A) - h(x-{Ay, A) = s'^ therefore defines a unique mapping z = g(s) with g'(0) = 
•^2/My(x_(A)) > 0 that is conformal on a complex neighborhood of i = 0 containing a square 
max{|Re{i)|, |Im{i)|) < 6' whose conformal image in turn contains a neighborhood of the form 
|z - x_(A)| < (5 for some 6 - 6{A) > 0 (independent of e). If -L < x < x_(A) - <5, we integrate 
by parts using (A^ + to obtain 




V27re 


£ 


-ihc(z)/c 


dz 




(C.ll) 


for some K' - K'(d) independent of both e and x 6 [-L, x_(A) - 5]. If instead x_(A) + d < x < 
x(A), we first integrate from -L to x(A) and then back to x, and obtain by the same argument 


ap,e(x) - ap,A-L) ■ 


1 


V27re 


I 


■x(A) 


1 


yj2ne 


I 


x(A) 


(C.12) 


The integral on the left-hand side (over -L < z < x(A)) is independent of x, and by the method 
of stationary phase (appropriately generalized as described in Remark [5T| using Id^ - A| < Ce) 
it is Oi -sje) due to the simple critical point at z = X-(A). Therefore, combining ( |C.l l| l-( |C.12| i 
shows that |a' ^(x)| < K also holds for -L < x < X-(A) - 6 and for x_(A) -Hd < x < x(A). Finally, 
for |x - x_(A)| < 5 we write 


ap,Ax) - ap,A-L) = p— 

ylne 


1 dz 1 


X 


z - z, 


where 




V27re 

g{s) - Zp 


rg W 

l’^kp,,{s)ds (C.13) 

Jo 


(C.14) 


The first term on the right-hand side of ( |C.13| l is independent of x and is (9(1) as e ^ 0 by 
(generalized) stationary phase. To analyze the second term, we assume that x > x_(A) (the case 
X < x_(A) is similar) which implies that g '(■*) > 0, and apply Cauchy’s Theorem in the square 
max{|Re{s)|, |Im{s)|) < S', on which kp^^{s) is analytic and bounded uniformly for small ebecause 
|df - A| < Ce, to deform the real path to a diagonal and a vertical segment; 


rg ‘(x) 

Jo ' 


e ‘^^^"kp^^is) ds - 


But because |kp,e(i)| < K", 

-{\-i)g-Hx) 


X 


X 


e "‘^^^kp^^{s)ds 


<\-i)g-'U) 


e ''^^^kp^fis) ds ■ 


r'"’ 

J(l-0s“'(x) 


e’'’’^l^kp^,{s)ds. (C.15) 
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and similarly 


rg (-t) , p 

e-'^'/^kp,,{s)ds <K" 
J(l-i)e-'(x) Jo 


^g-'M 


-~g-'(x)t/€ 


dt 


< K' 


rg 'W 

Jo ' 


~t^le 


ds < 


K” 


(C.17) 


Using ( |C.15| )-( |C.17[ ) in ( |C.13| l then shows that |fl^^(x)| < K holds also for \x - x_(A)| < 6. 

Property 3. Finally, the pointwise limit ( |60| follows from an application of the method of sta¬ 
tionary phase (or steepest descent) to J,x) (for x < x(A), x + x_(A)) and ^(x) (for x > x(A), 
X + x+(A)), taking into account Remark 5.1 and the fact that \A^ - A| < Ce. The integrals in 
( |C.7| l involve contributions from stationary phase points only if X-(A) < x < x+(A), and hence 
aj,^{x) 0 pointwise for x < x_(A) while a^Jx) —> 0 pointwise for x > x+(A), implying the 

same trivial limit for 4(x). For x e (x_(A), x+(A)), the integrals a^pjx) have a unique simple 
stationary phase point at z = x_(A), with the result that 




g-i;r/4g-!/!(x_(A);ie)/e 


+ 0 ( Vi), 


0 . 


(x-(A) - Zp) ^ju'^ix-iK)) 

Using this result in (|C.6|l and using Lemma C. 1 then gives, for x_(A) < x < x+(A), 


(C.18) 


lim4(x) = - 

5-^0 Am[)(x-(A)) 


z 

p=l 




x-(A) ■ 


m'(x-(A)) 


u'(x-(A)) 


z 


nlL/(zp - 


pA A-(^) - Zp) Wii.p{zp - Zi) 




c (x-(A)-z) UUJz-zi) 


■ dz 


(C.19) 


where C is a positively-oriented contour enclosing zi, ■. ■,Zp but excluding x-(A). Evaluating 
this integral by taking an exterior residue at z = x-(A) (no residue at oo) and comparing with 
( [T8] i and ( [T9| l establishes the limit ( |60| ) when x_(A) < x < x+(A). 
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